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LONG  TERM  GOALS 

Time  and  space  dependent  radiance  distributions  at  the  sea  surface  are  a  function  of  the  shape  of  the 
incident  distribution  on  the  surface,  modification  by  the  sea  surface  itself  from  topography  and 
transmission  characteristics,  and  alteration  by  the  Inherent  Optical  Properties  (lOPs)  of  the  surface 
ocean.  Our  long  term  goal  is  understanding  this  last  controlling  factor.  With  a  knowledge  of  the  lOPs, 
radiance  fields  can  be  directly  computed  from  the  incident  field  using  the  equation  of  radiative  transfer, 
now  embedded  in  commercially  available  code  (e.g..  Hydrolight). 

With  the  state  of  current  technology  and  methodologies,  the  primary  obstacles  in  understanding 
subsurface  lOPs  and  their  high-frequency  dynamics  are  a  lack  of  1)  volume  scattering  instrumentation, 
2)  comprehensive  inversion  models  linking  the  lOPs  with  the  ambient  particle  fields  including  bubbles 
(models  which  in  many  cases  will  require  input  dependent  on  1),  and  3)  suitably  stable,  non-intrusive 
platforms  to  sample  the  subsurface  ocean.  The  first  two  challenges  are  addressed  in  this  project. 

OBJECTIVES 

There  are  two  broad  objectives  for  this  project: 

1)  To  develop  an  in-situ  volume  scattering  function  device  measuring  volume  scattering  from  10°  to 
170°  at  10°  intervals  and  sampling  rates  of  1  s  ’  or  better  to  sample  the  VSF  in  near-surface  waters;  and 

2)  To  develop  and  refine  lOP  inversion  models  to  resolve  particle  field  characteristics  on  small 
spatial  (cm’s)  and  temporal  (<1  s)  scales  in  near-surface  waters. 
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APPROACH 


Our  design  for  the  VSF  device  is  illustrated  in  Figure  1.  The  device  is  called  MASCOT  (Multi- Angle 
scattering  Optical  Tool).  The  source  beam  is  a  30  mW  658  nm  laser  diode  expanded  with  a  Gallilean 
2X  beam  expander  to  an  approximately  3  mm  X  8  mm  elliptical  shape.  A  wedge  depolarizer  is  used  to 
provide  the  unpolarized  light  needed  for  VSF  determinations.  Seventeen  independent  silicon  diode 
detectors  spaced  in  a  semicircle  10  cm  around  the  sample  volume  measure  volume  scattering  from  10° 
to  170°  at  10°  intervals.  The  total  pathlength  for  all  scattering  measurements  is  20  cm.  Independent 
detectors  allow  resolution  of  the  VSF  without  any  moving  parts  and  time-consuming  scanning. 
Additionally,  each  detector  can  be  optimized  for  its  specific  dynamic  range.  Detector  field-of- views 
(FOVs)  range  from  0.8°  to  5°  for  the  different  detectors,  with  the  narrowest  FOVs  associated  with  the 
detectors  measuring  scattering  in  the  forward  direction.  Using  proprietary  electronics,  a  20  Hz 
sampling  rate  for  all  channels  has  been  achieved  while  maintaining  a  worst  case  signahnoise  of  300: 1 . 
Relatively  fast  sampling  rates  are  important  in  resolving  VSFs  in  the  highly  dynamic  ocean  subsurface. 

Polarized  VSFs  have  now  been  successfully  collected  with  the  addition  of  filter  mount  placed  in  front 
of  the  source  beam.  A  linear  polarizer  is  used  to  obtain  scattering  from  a  vertically  and  horizontally 
polarized  source  (in  terms  of  the  Mueller  scattering  matrix  elements,  (SI  l+S12)/2  and  (SI  l-S12)/2  are 
measured,  so  that  SI  1  and  S12  may  be  derived).  Adding  polarized  scattering  increases  the  amount  of 
information  on  particle  characteristics  we  are  collecting,  and  is  expected  to  improve  our  ability  to 
discriminate  different  particle  types  (both  the  number  of  subpopulations  and  the  accuracy  of  individual 
determinations).  The  degree  of  linear  polarization  (S12/S1 1)  is  most  dependent  on  the  degree  of 
sphericity,  particle  size,  and  refractive  index  composition. 

For  VSF  inversion  modeling,  we  are  extending  the  capabilities  of  existing  models  (Twardowski  et  al. 
2001;  Twardowski  and  Zaneveld,  2004;  Zhang  et  al.  2005)  by  incorporating  input  from  new  VSF 
measurements  and  by  adding  bubble  particle  populations  (clean  and  coated)  in  the  models.  Candidate 
phase  functions  for  particle  subpopulations  are  fit  to  measured  VSFs  using  a  least-squares 
minimization  procedure.  These  phase  functions  can  be  obtained  theoretically  using  Mie  theory,  DDA, 
or  FDTD  techniques,  or  experimentally  in  controlled  laboratory  conditions. 

The  MASCOT  has  now  been  extensively  deployed  concurrently  with  the  commercially  available  near¬ 
forward  VSF  device  LISST  (Sequoia  Inc.)  in  order  to  capture  the  VSF  with  good  resolution  from  ~0.1 
degrees  to  170  degrees.  Our  deployment  package  also  contains  CTD,  AC9  (or  ACS),  and  various  ECO 
scattering  sensors,  all  integrated  in  a  real-time  vertical  profiling  system  (Fig.  2). 

Deployments  for  RaDyO  have  taken  place  in  collaboration  with  other  RaDyO  investigators  off  Scripps 
Pier  in  January  2008,  in  Santa  Barbara  Channel  in  September  2008,  and  are  ongoing  off  Hawaii  at  this 
time.  For  these  RaDyO  deployments,  we  have  also  intermittently  integrated  two  additional  sensors:  1) 
a  high  sampling  rate  fish-eye  lens  radiometer  currently  being  developed  by  Marlon  Lewis  and  Scott 
McLean  of  Satlantic,  and  2)  a  bubble  acoustic  resonator  developed  by  Svein  Vagle,  David  Farmer,  and 
Helen  Czerski.  For  the  Scripps  Pier  experiment,  a  mobile  deployment  platform  was  configured  so  that 
measurements  could  be  made  along  the  length  of  the  pier,  with  the  expectation  that  higher 
concentrations  of  bubbles  would  be  observed  nearing  the  surfzone.  For  subsequent  field  work  in  SBC 
and  Hawaii,  we  are  deploying  our  sensor  package  off  the  RA/  Kilo  Moana. 
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Figure  1.  Oblique  view  illustration  of  the  MASCOT.  The  VSF  is  resolved  from  10  to  170  degrees  in 
10  degree  intervals.  Detectors  are  wedge  shaped  and  arranged  in  a  semi-circle  on  an  aluminum 
frame  to  minimize  reflections  and  perturbation  of  the  water  sample  in  the  remote  volume  (center  of 
semi-circle).  The  source  assembly  includes  a  30  mW 658  nm  laser  diode,  reference  detector,  beam 
expander,  and  wedge  depolarizer.  Wiring  from  all  the  detector  modules  and  the  source  module 

feeds  to  a  data  handling  unit. 


WORK  COMPLETED 

•  An  lOP  package  including  the  MASCOT  and  LISST  was  deployed  off  Scripps  Pier,  in  the  Santa 
Barbara  Channel,  and  off  Hawaii  during  collaborative  RADYO  exercises. 

•  The  MASCOT  was  adapted  to  additionally  measure  linear  polarized  scattering  elements  to  derive 
degree  of  linear  polarization  (DOLP)  as  a  function  of  angle. 

•  An  acoustic  bubble  resonator  (Vagle  and  Farmer)  and  a  fish-eye  lens  camera  (Lewis)  were  successfully 
integrated  on  the  MASCOT  lOP  package  and  deployed  periodically  for  all  exercises. 

•  The  MASCOT  lOP  package  was  also  tested  off  the  New  York  bight  in  May  and  July  in  association 
with  other  funded  work. 

•  A  revised  calibration  protocol  was  developed  and  implemented,  using  clean  techniques  and  accurately 
accounting  for  pure  water  absorption  along  the  optical  path. 

•  A  revised  calibration  and  correction  protocol  was  developed  for  the  LISST  to  account  for  ambient  light 
contamination,  turbulent  scattering,  and  vignetting 

•  A  second  MASCOT  prototype  with  a  larger  dynamic  range  suitable  for  high  turbidity  environments 
(such  as  the  surf  zone)  was  tested  in  the  lab  and  deployed  in  the  NY  bight. 

•  Inversion  of  MASCOT  VSFs  collected  off  Scripps  Pier  and  the  SBC  into  component  VSFs  associated 
with  subpopulations  of  the  bulk  particle  field  including  bubbles  has  been  refined  using  a  least-squares 
minimization  model  based  on  Mie  theory,  including  derivation  of  particle  concentration  and  size 
distribution  for  each  subpopulation.  We  are  currently  adapting  this  model  to  also  invert  for 
nonspherical  asymmetric  polyhedron  particles  more  analogous  to  natural  suspended  mineral 
populations. 

•  3  papers  are  in  prep  and  1  was  just  accepted  that  are  based  on  our  results  from  this  project. 
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RESULTS 


Fig.  2  shows  the  MASCOT  device  in  a  custom  cage  with  additional  optical  sensors  ready  for 
deployment.  Orientation  in  the  horizontal  plane  minimizes  any  shear  perturbation  of  the  sampled 
water  parcel  during  upcasts. 

Scripps  Pier.  Fig.  3  shows  scattering  results  from  our  deployment  of  the  MASCOT  off  Scripps  Pier  in 
January.  Typically,  waves  breaking  shoreward  of  the  sensor  package’s  location  produced  suspended 
sediment  and  bubble  plumes  that  became  entrained  in  a  seaward  rip  current  running  along  the  south 
side  of  the  pier  where  the  sensor  package  was  suspended.  Bubble  plumes  had  a  visible  reflectance 
expression  for  several  seconds,  while  discrete  suspended  sediment  plumes  were  clearly  observed 
lasting  for  up  to  minutes.  At  high  tide  (the  condition  at  the  beginning  of  the  Fig.  3  time  series),  with 
the  surf  zone  far  shoreward  of  our  position  (>100  m),  plumes  with  injected  air  bubbles  typically  did  not 
persist  long  enough  to  reach  our  sensors,  although  episodic  local  breaking  with  localized  bubble 
injection  did  occasionally  occur.  In  these  conditions,  periodic  plumes  of  suspended  sediment  entrained 
in  the  seaward  rip  current  were  clearly  visible  as  turbidity  maxima  with  a  distinct  VSF  character 
(greater  relative  backscattering  than  the  background).  As  the  tide  ebbed,  the  surf  zone  passed  directly 
through  our  sampling  region  so  that  plumes  of  both  bubbles  and  sediments  were  sampled.  Plumes 
dominated  with  bubble  particles  exhibited  a  well-known  enhancement  in  scattering  at  mid-angles,  from 
60  to  90  degrees  (e.g.,  Zhang  et  al.  2004).  This  is  the  first  time  this  mid-angle  enhancement  in 
scattering  due  bubbles  has  been  observed  in-situ.  The  20  Hz  sampling  rate  for  all  VSF  channels  allows 
sufficient  resolution  of  the  bubble  plumes  as  they  rapidly  evolve  in  time  and  space.  A  very  low 
background  scattering  condition  (particle  attenuation  less  than  0.3  m''  at  532  nm)  enhanced  our 
dynamic  range  in  resolving  intense  scattering  associated  with  episodic  bubble  and  sediment  plume 
generation. 

A  library  of  VSFs  generated  from  a  coated  Mie  theory  model  representing  180  initial  particle 
subpopulations  were  fitted  to  measured  VSFs  using  a  least-squares  minimization  technique  (see  Fig.  4 
for  details).  This  work  has  been  carried  out  in  collaboration  with  Xiaodong  Zhang  at  the  U.  of  N. 
Dakota.  Fig.  5  shows  inversion  results  for  combined  bubble  subpopulations  and  combined  mineral 
populations.  Results  are  encouraging  and  are  consistent  with  anecdotal  evidence  (what  we  saw)  with 
respect  to  patches  of  sediment  and/or  bubbles.  Fig.  6  shows  the  inversion  results  for  the  2  primary 
bubble  subpopulations,  and  Fig.  7  shows  results  for  the  primary  mineral  subpopulations.  Fig.  8  shows 
acoustics  data  collected  with  the  Svein/Farmer  bubble  resonator  concurrently  with  the  MASCOT  VSF 
data  for  the  last  40  minutes  of  the  Fig.  3  time  series.  Good  agreement  between  patterns  in  optical 
scattering  from  bubbles  (positive  displacements  of  the  P(60)/|3(120)  relationship  relative  to  the 
background)  and  acoustic  attenuation  at  multiple  frequencies  was  observed,  but  the  agreement  was 
even  better  when  comparing  with  the  large  bubble  subpopulation  from  the  inversion  results.  Fig.  9 
shows  the  size  distributions  from  the  inversion  results,  with  an  aggregate  distribution  consistent  with 
the  familiar  Junge-type  distribution  in  the  optically  relevant  size  range  (~0.1  to  ~40  pm).  Stable  model 
solutions  (only  possible  to  assess  with  a  very  large  number  of  VSFs)  consistent  with  anecdotal 
observations  of  particle  composition,  coupled  with  validation  based  on  acoustic  measurements  of 
bubble  distributions  and  fractional  and  aggregate  size  distributions  agreeing  with  theory  and 
observations,  provides  strong  evidence  that  our  inversion  results  are  environmentally  meaningful. 

Santa  Barbara  Channel.  During  the  recent  Santa  Barbara  Channel  experiment,  over  150  casts  or  time 
series  were  collected,  including,  for  the  first  time,  polarized  scattering  elements  with  quasi-concurrent 
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measurements  of  dark  offsets  to  ensure  the  greatest  accuracy  possible.  Plots  of  attenuation, 
chlorophyll,  backscattering  ratio,  and  bulk  refractive  index  for  the  entire  cruise  are  provided  in  Fig.  10- 
13.  Final  corrected  VSF  data  from  MASCOT  and  LISST  sensors  are  provided  in  Fig.  14  for  sample 
casts.  Optical  scattering  from  seawater  refractive  index  discontinuities  is  always  present  in  the  LISST 
near-forward  measurements  (<~0.1  deg).  Vertical  profiles  of  MASCOT  scattering  are  shown  in  Fig. 

15.  A  comparison  of  total  scattering  derived  from  VSF  measurements  and  total  scattering  from  the  ac9 
showed  good  agreement  (Fig.  16).  Degree  of  linear  polarization  (DOLP)  results  for  2  casts  are  shown 
in  Fig.  17  and  are  consistent  with  the  previous  bench  top  results  of  Voss  and  Fry  (1984).  The 
maximum  angle  of  DOLP  was  typically  around  94  deg.  (Fig.  18).  Mie  theory  results  modeling  DOLP 
for  bubble  populations  showed  a  markedly  different  function  than  is  typical  for  particles  (Fig.  19),  with 
minimal  DOLP  (i.e.,  ~0)  between  0  and  80  deg.  Such  a  unique  DOLP  signature  could  serve  as  an 
effective  marker  for  bubbles  in-situ.  This  was  testing  in  SBC  with  a  set  of  subsurface  time  series  under 
conditions  of  actively  breaking  waves  (Fig.  20).  When  healthy  bubble  populations  were  present  in 
both  horizontal  and  vertically  polarized  measurements,  results  were  in  reasonable  agreement  with 
theory  (Fig.  21).  Another  example  from  the  NY  Bight  with  white-capping  conditions  also  shows 
reasonable  agreement  with  theory  (Fig.  22).  These  data  underscore  the  need  to  switch  between 
polarization  states  rapidly  in-situ  (as  opposed  to  consecutive  time  series  of  different  scattering 
polarization  states),  which  we  are  currently  working  on.  We  are  also  working  on  developing  compact, 
low-power  ECO-type  sensors  for  automated  platforms  to  make  VSF  measurements  suitable  for 
discriminating  the  mid-angle  enhancement  due  to  bubbles  (Fig.  23). 

Hawaii.  The  shape  of  the  mean  MASCOT  VSF  for  the  entire  cruise  was  typical  of  previous  data  (Fig. 
24;  also  see  Figs.  27,  28).  Largest  standard  deviations  were  located  at  mid-angles  due  to  the 
intermittent  presence  of  bubbles  at  the  surface.  Mean  attenuation  and  absorption  spectra  are  plotted  in 
Fig.  25  and  are  typical  for  the  region.  The  turbulent  component  of  LISST  scattering  matched  refractive 
index  discontinuities  in  the  water  column  (Fig.  26),  consistent  with  previously  reported  “schleiren” 
phenomena  (Mikkelsen  et  al.  2006). 

VSF  shape  in  the  backward  direction  was  remarkably  consistent  among  the  several  million  MASCOT 
VSFs  collected  in  8  field  sites  exhibiting  a  wide  range  in  water  types  (Figs  27,  28). 

IMPACT/APPLICATIONS 

Progress  and  results  represent  important  steps  toward  the  development  and  vetting  of  a  multi-angle,  in¬ 
water  VSF  device.  Knowledge  of  the  Inherent  Optical  Properties  including  the  VSF  can  be  used  to 
predict  and  optimize  the  performance  of  a  host  of  Naval  operations  that  rely  on  divers,  cameras,  laser 
imaging  systems,  and  active  and  passive  remote  sensing  systems.  These  include  mine 
countermeasures,  harbor  security  operations,  debris  field  mapping,  anti-submarine  warfare,  and  search 
and  salvage  operations. 
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Figure  2.  MASCOT  VSF  device  mounted  with  other  optical  sensors  in  a  custom  cage 
designed  to  sample  the  subsurface  domain. 
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Figure  3.  A  90  minute  time  series  of  the  volume  scattering  function  resolved  at  60  and  120  degrees  with  the 
MASCOT  off  Scripps  Pier.  At  the  beginning  of  data  collection  the  surf  zone  was  shoreward  of  the  sensor 
package  and  beginning  to  ebb.  Breaking  waves  generated  mixed  bubble  and  sediment  plumes  which 
became  entrained  in  a  seaward  rip  current  along  the  pier.  By  the  time  these  plumes  reached  the  sensor 
package  (tens  of  seconds),  the  visible  bubble  influence  had  dissipated,  leaving  the  strong  sediment 
contribution.  Maximum  wave  sets  broke  onshore  every  ~18  minutes.  Scattering  plumes  from  the  first  4  sets 
were  dominated  by  sediment,  exhibiting  a  decreased  P(60)/P(120)  relationship  due  to  the  sediment 
dominated  plumes  having  a  higher  refractive  index  than  the  background  particle  population  (e.g., 
Twardowski  et  al.  2001).  Conversely,  when  the  tide  had  sufficiently  ebbed  so  that  the  surf  zone  was 
adjacent  to  the  sensor  package  and  active  bubble  plumes  were  being  directly  observed,  scattering  plumes 
exhibited  an  increased  P(60)/p(120)  relationship  due  to  the  enhancement  of  mid-angle  scattering  due  to 

bubbles  (e.g.,  Zhang  et  al.  2004). 
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Figure  4.  Phase  functions  computed  from  Mie  theory  used  in  the  VSF  inversion  algorithm.  From  180 
initial  phase  functions  (particle  parameters  included  5  bulk  refractive  indices,  12  modal  diameters  for  a 
lognormal  distribution,  3  standard  deviations  for  the  distributions,  and  a  bubble  coating  thickness  scaled  at 
2%),  59  were  determined  to  be  >10%  rmse  different.  After  binning  the  MASCOT  VSF  time  series  data  from 
Scripps  to  1  Hz  (N  =  5395),  19  phase  functions  (plotted  here)  from  the  library  were  picked  by  the  least- 

squares  fitting  algorithm  to  fit  these  measured  VSFs. 
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Figure  5.  VSF  inversion  results  for  the  Scripps  Pier  time  series j  including  panels  of  P(60)/  P(120)  (see  Fig. 
3)y  the  contribution  from  all  bubble  subpopulations  to  total  scatterings  and  the  contribution  from  all 
mineral  populations  to  total  scattering.  The  5  primary  wave  sets  are  marked  in  gray.  Results  show  that 
increases  in  scattering  during  the  first  4  wave  sets  were  primarily  due  to  suspended  sediments.  After  the  4^ 
wave  sets  increases  in  scattering  were  due  to  bubbles  and  sediments. 
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Figure  6.  VSF  inversion  results  for  the  two  dominant  bubble  subpopulations,  a  small  subpopulation 
with  modal  diameter  and  standard  deviation  of  0.2  pm,  and  a  larger  subpopulation  with  modal  diameter 
and  standard  deviation  of  2  pm  and  0.6  pm,  respectively.  Phase  functions  of  the  two  subpopulations  are 
shown  in  the  inset  of  each  panel.  The  mid-angle  enhancement  due  to  larger  bubbles  is  evident  for  the 
larger  bubble  population.  There  is  evidence  of  the  small  bubble  subpopulation  in  the  inversions  from 
the  particle  maxima  observed  during  the  first  4  wave  sets.  This  appears  to  be  evidence  of  controversial 
persistent  small  bubbles  (<10  pm)  in  these  “aged”  plumes  of particles  (order  of  minutes  old).  Later  in 
the  time  series  when  the  waves  are  breaking  close  to  the  instrumentation,  frequent  and  substantial 
peaks  in  scattering  from  both  small  and  large  bubbles  are  observed.  Some  succession  can  be  observed 
between  the  large  bubble  population  and  the  more  persistent  small  bubble  population  on  the  order  of 

seconds. 
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Figure  7.  VSF  inversion  results  for  the  dominant  sediment  subpopulations,  plotted  in  5  subpanels, 
ordered  from  smallest  to  largest.  Results  show  the  most  important  sediment  subpopulations  are  the 

smallest  two. 
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Figure  8.  Inversion  results  from  the  large  bubble  population  (upper  panel)  plotted  with  relative 
acoustic  attenuation  from  a  bubble  resonator  (lower  panel).  Patterns  were  very  consistent,  more 
consistent  than  observed  when  comparing  the  acoustic  attenuation  with  positive  responses  in 
P(60)/p(120)  from  the  MASCOT  (also  plotted  in  lower  panel).  The  response  in  the  simple  ratio 
P(60)/p(120)  can  get  complex  when  both  bubbles  and  sediment  predominate.  As  a  result,  on  some 
occasions  this  ratio  disagrees  with  the  acoustic  attenuation  (e.g.,  at  minute  66.5),  whereas  the  more 
rigorous  VSF  inversion  results  show  better  agreement.  The  time  series  shown  here  is  a  subset  of  the 
data  shown  in  Fig.  3,  including  only  wave  sets  4  and  5.  Note  that  the  turbidity  maximum  observed 
around  59-60  minutes  in  Fig.  3  is  associated  with  a  decreasing  P(60)/p(120),  no  acoustic  response,  and 
only  a  small  response  in  the  large  bubble  inversion  results.  This  is  consistent  with  the  previous 
interpretation  that  this  turbidity  plume  is  dominated  by  sediments.  A  particularly  close  agreement  was 
observed  between  the  inversion  results  and  acoustic  attenuation  at  the  frequency  specific  for  143um 
diameter  bubbles.  Note  that  even  though  the  two  sensors  were  on  the  same  instrument  package,  sample 

volumes  were  separated  by  about  60  to  70  cm. 
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Figure  9.  Size  distributions  of  dominant  subpopulations  with  relative  refractive  indices  of  0. 75 
(bubbles),  1.16  (quartz-like  minerals),  and  1.22  (clay-like  minerals).  The  aggregate  size  distribution  (in 
black)  is  modeled  well  with  a  hyperbolic  differential  Junge-type  model  with  slope  3. 7  in  the  range  that 
can  be  resolved  by  sizing  devices  such  as  a  coulter  counter.  Many  previous  studies  looking  at  oceanic 
size  distributions  have  observed  that  the  Junge  distribution  is  a  good  model,  which  serves  as  an  element 

of  validation  for  our  inversion  results. 
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Figure  10.  Non-water  beam  attenuation  coefficient  at  650  nm,  cpg(650)y 
from  Santa  Barbara  Channel  2008. 
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Figure  11.  ac-9  derived  chlorophyll  from  Santa  Barbara  Channel  2008. 
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Figure  12.  Backscattering  ratio  at  657  nm,  bbp(657)/bp(657),  from  Santa  Barbara  Channel  2008. 
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Figure  13.  The  bulk  refractive  index  estimated  according  to  Twardowski  et  al.  (2001) 

from  Santa  Barbara  Channel  2008. 
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Figure  14.  Example  VSFs  plotted  as  a  function  of  angle  from  120  m  vertical  profiles  collected  with 
MASCOT  and  LISST  sensors  in  SBC.  Agreement  between  MASCOT  and  LISST  sensors  in  the  region 
of  overlap  (around  10  deg)  is  excellent.  Sensor  calibrations  are  completely  independent.  Optical 
scattering  from  turbulent  refractive  index  discontinuities  is  present  in  all  near-forward  LISST  data  (less 
than  ~0.1  deg).  In  some  casts  such  as  SBC  149,  a  mid-angle  enhancement  due  to  bubbles  can  be 

observed  in  some  (subsurface)  VSFs. 
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Figure  15.  MASCOT  VSF  data  plotted  as  function  of  depth  for  all  1 7  angles  except  30  deg  (faulty 
detector).  A  surface  particle  maximum  was  always  present.  A  deep  maximum  of  resuspended  sediments 
was  sometimes  present.  The  clearest  water  was  observed  around  80  m. 
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Figure  16.  Good  agreement  was  observed  between  total  particle  scattering  derived  from  integrating  the 
full  volume  scattering  function  using  the  acceptance  angle  of  the  ac9  (0.93  deg)  as  the  minimum  angle 
and  total  particle  scattering  measured  with  an  ac9.  These  data  are  from  one  120  m  cast  on  Sept,  day 

19.69. 
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Figure  1 7.  MASCOT  degree  of  linear  polarization  (DOLP)  as  a  function  of  angle  for  2  example 
profiles,  the  second  of  which  sampled  a  deep  water  particle  maximum  comprised  of  resuspended 
sediments  (with  corresponding  minimum  in  DOLP).  Shapes  are  generally  consistent  with  Voss  and  Fry 
(1984),  where  maxima  typically  occur  around  0.8  and  at  angles  just  larger  than  90  deg;  near  0  and  180 

DOLP  converges  on  zero. 
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Figure  18.  Angles  for  MASCOT  DOLP  maxima  for  the  profiles  from  Fig.  1 7,  estimated  with  a  cubic 

spline  fitting  algorithm. 
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Figure  19.  Theoretically  derived  DOLP  for  lognormally-distributed  bubble  populations.  The  mean 
radius  varied  according  to  ln(l:  1:100)  and  the  standard  deviation  for  the  distribution  varied  according 
to  10%  of  the  mean  radius.  Size  distributions  are  plotted  in  inset  graph.  Results  show  that  DOLP  for 
bubble  populations  can  be  expected  to  be  very  different  in  shape  than  DOLP  for  oceanic  particles  (e.g., 

see  Fig  1 7). 
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Figure  20.  Subsurface  time  series  of  MASCOT  VSF  scattering  at  80  deg.  (location  of  mid-angle 
enhancement  due  to  bubbles)  for  an  unpolarized  source  (top  panel),  a  vertically  polarized  source 
(middle  panel),  and  a  horizontally  polarized  source  (bottom  panel).  Time  series  were  collected 
consecutively,  requiring  retrieval  of  the  instrument  package  and  manual  polarization  switching  between 
each  series.  Data  was  collected  during  20-25  knt  winds  accompanied  by  extensive  white  capping.  All 
observed  enhancement  above  the  baseline  in  each  time  series  is  due  to  injected  bubble  particle  fields. 
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Figure  21.  DOLP  computed  from  all  combinations  of  1  minute  time  averages  from  the  time  series 
plotted  in  Fig.  20.  In  cases  where  healthy  bubbles  are  observed  in  both  the  vertically  and  horizontally 
polarized  source  measurements,  the  resulting  DOLP  showed  reasonable  agreement  with  the  theoretical 
expectation.  For  those  cases  where  bubbles  were  in  low  abundances,  the  resulting  DOLP  showed  a 
more  typical  shape  consistent  with  background  particles  (e.g.,  see  Fig.  1 7).  For  those  cases  where 
bubbles  were  moderate  abundance  relative  to  background  particles,  a  DOLP  in  between  these  DOLPs 
was  observed.  However,  when  healthy  bubbles  were  observed  in  one  polarization  measurement  and  not 
the  other,  unrealistic,  strongly  negative  DOLPs  were  observed. 
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Figure  22.  DOLP  results  observed  in  NY  Bight  under  extensive  whitecapping.  Again,  reasonable 
agreement  between  measurement  and  theoretical  expectation  is  observed. 
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Figure  23.  Newly  developed  low  power,  compact  ECO-type  sensor  for  discriminating  the  mid-angle 
scattering  enhancement  characteristic  of  bubbles  greater  than  a  several  /jm  in  diameter.  Sensor  is  suitable 

for  deployment  on  compact  autonomous  platforms. 
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Figure  24.  Mean  and  standard  deviation  of  MASCOT  non-seawater  VSF  for  Hawaii  cruise  AUG- 
SEP  2009.  Note  largest  standard  deviations  are  associated  with  middle  angles  where  intermittent 

bubbles  exhibit  enhanced  scattering. 
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Figure  25.  Mean  non-seawater  attenuation  (Cpg)  and  absorption  (apg) 
spectra  for  Hawaii  cruise  AUG-SEP  2009. 
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Figure  26.  Seawater  refractive  index  gradient  as  a  function  of  depth  and  optical  turbulent  scattering  in  the 
near-forward  measured  with  the  LISST.  A  strong  gradient  in  refractive  index  is  observed  around  60  m  with 
a  concomitant  increase  in  turbulent  scattering  derived  from  the  LISST  data. 
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Figure  27.  Mean  and  standard  deviation  of  MASCOT  phase  function  data  in  the  backward  direction  for 
several  million  VSFs  collected  during  field  exercises  in  8  sites  ranging  widely  in  water  type  from  clear 
waters  off  Hawaii  and  in  the  Ligurian  Sea  to  turbid  waters  of  NY  Harbor  (circles  with  dashed  line  in  upper 
panel).  The  coastal,  clear,  and  turbid  water  phase  function  data  of  Petzold  (1972)  are  also  plotted  in  the 
upper  panel,  with  the  turbid  measurement  showing  the  best  agreement.  The  smallest  standard  deviation  in 
phase  function  was  observed  at  110  and  120  deg.  Fournier-Forand  phase  functions  exhibit  a  similar  shape 
to  our  averaged  phase  function  in  the  backward  and  also  show  smallest  deviation  between  110  and  120  deg. 
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Figure  28.  The  factor  x  is  used  to  convert  a  VSF  measurement,  P,  to  a  backscattering  coefficient 
according  to  bbp  =  Xn  271  PjfO),  and  is  most  effective  at  the  angle  where  the  backward  phase  function  shows 
the  smallest  variability  in  an  extensive  data  set.  The  xus  a  function  of  angle  is  plotted  for  our  data  set  in 
the  upper  panel  along  with  previous  results  from  isolated  regions  collected  with  a  Ukrainian  device  called 
the  MVSM  (Boss  and  Pegau  2001;  Berthon  et  al.  2007;  and  Chami  et  al.  2005).  Our  results  agree 
reasonably  well  with  previous  results  except  for  Chami  et  al.  MVSM  results  in  general  exhibit  lower  x 
values  in  the  near-backward,  perhaps  consistent  with  a  reflection  problem,  which  is  always  a  challenge  for 
these  types  of  measurements.  The  lower  panel  shows  the  %variability  for  our  phase  function  data  along 
with  the  MASCOT  electronic  noise  estimated  from  in-situ  dark  count  measurements.  Remarkably,  the 
baseline  noise  is  comparable  to  the  ~2%  variability  observed  the  phase  function  around  110  and  120 
degrees,  indicating  that  almost  all  the  variability  observed  is  instrument  noise  and  the  subsequent  “real” 
variability  in  the  phase  function  at  these  angles  is  practically  negligible. 
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TRANSITIONS 


We  expect  that  our  efforts  in  developing  an  in-water  VSF  device  and  associated  inversion  techniques 
to  better  understand  particle  dynamics  in  natural  waters  will  lead  to  transition  as  operational  tools  for 
the  fleet  and  the  oceanographic  research  community  in  the  future. 

RELATED  PROJECTS 

This  effort  is  related  to  several  ongoing  efforts  by  the  PI  to  develop  optical  sensors  and  associated 
biogeochemical  inversion  techniques  to  improve  our  understanding  of  the  oceanic  environment. 
Current  ongoing  related  projects  include: 

•  investigating  the  dynamics  of  scattering  by  subsurface  bubble  populations  and  other  particles  in 
the  S.  Ocean  (NASA,  Twardowski;  project  lead  PI  H.  Dierssen); 

•  investigating  the  underlying  controls  of  biological  camouflage  responses  in  dynamic 
underwater  optical  environments  (MURI  collaboration,  Twardowski;  URI  Pis  J.  Sullivan  and 
B.  Seibel;  project  lead  PI  M.  Cummings); 

•  developing  improved  remote  sensing  water  quality  algorithms  for  coastal  waters  (NASA, 
Twardowski;  project  lead  PI  Z-P.  Lee); 

•  developing  compact,  low  power  sensing  tools  for  ocean  observing  platforms  (ONR  SBIR, 
Twardowski); 

•  developing  a  microscopic  holographic  camera  for  optically  relevant  particles  (NOPP,  project 
lead  PI  J.  Sullivan); 

•  developing  optical  prediction  models  for  the  surfzone  (CEROS,  Twardowski;  project  lead  PI  G. 
Chang);  and 

•  developing  a  surfzone  drifter  measuring  optical  attenuation  and  scattering  (ONR  SBIR, 
Twardowski). 
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